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Abstract 

We present a method how to estimate from experimental data of a tur¬ 
bulent velocity field the drift and the diffusion coefficient of a Fokker-Planck 
equation. It is shown that solutions of this Fokker-Planck equation reproduce 
with high accuracy the statistics of velocity increments in the inertial range. 
Using solutions with different initial conditions at large scales we show that 
they converge. This can be interpreted as a signature of the universality of 
small scale turbulence in the limit of large inertial ranges. 


1 Introduction 

The common picture of fully developed local isotropic turbulence is that the ve¬ 
locity field is driven by external fields on large scales. By this driving energy is 
fed into the system at scales larger than the integral length Lq. A cascading pro¬ 
cess will transport this energy to smaller and smaller scales until at the viscous 
length scale r] the injected energy is finally dissipated by viscous effects [1, 2]. It 
is commonly believed that this picture of a cascade leads to the universality of the 
statistical laws of small scale turbulence. 

The standard quantity to study the statistics of turbulent fields is the longitu¬ 
dinal velocity increment of the length scale U defined as Vi = u(x + Li) — u(x), 
where u denotes the velocity component in the direction of the separation Li. x 
is a selected reference point. 

Two common ways have been established to analyze the statistical content of 
the velocity increments. On one side the structure functions have been evaluated 
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by supposing a scaling behavior 


< (vi) q >= Jdvi(vi) q P(vi,Li) ocL^, (1) 

where P(vi, Li), denotes the probability density function (pdf) for v , at the length 
scale Li. Using the recently proposed so-called extended selfsimilarity [3] it has 
become possible to evaluate the characterizing scaling exponents C, q quite accu¬ 
rately c.f. [4, 5]. On the other side one is interested to parameterize directly the 
evolution of the probability density functions (pdf) P(i>,:,L,;) c.f. [6, 7]. 

One major challenge of the research on turbulence is to understand small scale 
intermittency which is manifested in a changing form of P(vi,Li) or equivalently 
in a nonlinear g-dependence of the scaling exponents £ q , provided that the scaling 
assumption is valid. 

Although it is well known that, due to the statistic of a finite number of data 
points (let say 10 7 data), it is not possible to determine accurately scaling expo¬ 
nents (g for q > 6 [8], and that there are different experimental indications that 
no good scaling behavior is present [9, 10], for long time the main effort has been 
put into the understanding of the q-dependence of C, q (for actual reviews see [2, 5]. 
That there have been less attempts to analyze directly the pdfs may be based on 
the fact that up to now the scaling exponents are regarded as the simplest reduc¬ 
tion of the statistical content and that this analysis does not depend on model 
assumptions. In contrast to this, proposed parameterizations of the form of the 
pdfs [6, 7], although they are quite accurate, are still based on some additional 
assumptions on the underlying statistics. Based on the recent finding that the tur¬ 
bulent cascade obeys a Markov process in the variable L.j and that intermittency 
is due to multiplicative noise [10, 11, 12], we show in Section III that it is possible 
to estimate from the experimental data a Fokker-Planck equation, which describes 
the evolution of the pdfs with L*. We show that this Fokker-Planck equation 
reproduces accurately the experimental probability densities P(vi,Li) within the 
inertial range. Thus an analysis of experimental data is possible which quantifies 
the statistical process of the turbulent cascade, and which neither depends on scal¬ 
ing hypotheses nor on some fitting functions [13]. Having determined the correct 
Fokker-Planck equation for an experimental situation, in section VI we show how 
solutions of this equation with different large scale pdfs will converge to universal 
small scale statistics. This finding gives evidence for the universality of small scale 
turbulence. 


2 Experimental data 

The results presented here are based on 10 7 velocity data points. Local veloc¬ 
ity fluctuations were measured with a hot wire anemometer (Dantec Streamline 
90N10) and a hot wire probe (55P01) with a spatial resolution of about 1 mm. 
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The sampling frequency was 8 kHz. The stability of the jet was verified by mea¬ 
surements of the self-similar profiles of the mean velocity according to [14]. The 
turbulence measurements were performed by placing the probe on the axis of a 
free jet of dry air developing downwards in a closed chamber of the size 2m x lrn 
xlm. To prevent a disturbing counterflow of the outflowing air, an outlet was 
placed at the bottom of the chamber. The distance to the nozzle was 125 nozzle 
diameters D. As nozzle we used a convex inner profile [15] with an opening section 
of D = 8 mm and an area contraction ratio of 40. Together with a laminarizing 
prechamber we have achieved a highly laminar flow coming out of the nozzle. At a 
distance of 0.25 D from the nozzle no deviation from a rectangular velocity profile 
was found within the detector resolution. Based on a 12bit A/D converter reso¬ 
lution no fluctuations of the velocity could be found. The velocity at the nozzle 
was 45.5 m/s corresponding to a Reynolds number of 2.7 10 4 . At the distance of 
125 D we measured a mean velocity of 2.25m/s, a degree of turbulence of 0.17, an 
integral length Lq = 67mm, a Taylor length 0 = 6, 6mm (determined according to 
[16]), a Kolmogorov length r] = 0.25mm, and a Taylor Reynolds number Re = 190. 

The space dependence of the velocity increments Vi was obtained by the Taylor 
hypotheses of frozen turbulence. For the structure functions we found a tendency 
to scaling behavior for Lq > Li > 0. Intermittency clearly emerge as Li —► 0, as 
shown in Fig.l by the different form of the pdfs for different scales Li. 

3 Measurement of Kramers-Moyal coefficients 

Next, we show how to determine from experimental data appropriate statisti¬ 
cal equations to characterize the turbulent cascade. The basic quantity for this 
procedure is to evaluate the cascade by the statistical dependence of velocity 
increments of different length scales at the same location x [17]. Either two- 
increment probabilities p{v 2 , L 2 ,Vi,L]) or corresponding conditional probabilities 
p{v 2 ,L 2 \v 1 ,L 1 ) are evaluated from the whole data set. (Here we use the conven¬ 
tion that Tj+i < Li.) Investigating the corresponding three-increment statistics 
we could provide evidence that the evolution of these statistics with different Li 
fulfill the Chapman-Kolmogorov equation [18, 10, 11]. Furthermore evidences of 
the validity of the Markov process are shown as long as the step size between 
different Li is larger than an Markov length L mar , which is in the order of some p 
[12]. Thus we know [18] that the evolution of the increments Vi is given by a mas¬ 
ter equation without the involvement of memory functions and that the evolution 
of the pdfs P(ui, Li) are described by a Kramers-Moyal expansion: 

7 OO r\ 

- P(v i ,L i ) = Y / [-^] k D {k) (yi, L i) P(vi,Li). (2) 

dLi f—' out 
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Figure 1: Probability density functions P(v, L) determined from the experimental 
data (bold symbols) and calculated pdfs (lines) by the numerical iteration with the 
experimentally determined Fokkcr-Planck equation (5) and (6). The length scales 
A = ln(Lo/L) were: 0, 0.52, 1.04, 1.56 and 2.24 going from up to down. For the 
numerical iteration as initial condition an empirical fitting function for the large 
scale pdf was used which is shown by a broken line. For clarity of presentation the 
pdfs were shifted in y-direction. 


Where the Kramers-Moyal coefficients D^ k \vi, Li) are defined by the following 
conditional moments: 

M (k \v i,Ti,(5) = ^ Jdv 2 (v 2 ~ vi) k p{v 2 ,L 2 \v 1 ,L 1 ) 

D^ k \v u L r) = lim (iq, L\,6), (3) 

o —»0 

where S = L\ — L 2 . Knowing the Kramers-Moyal coefficients, and thus the evo¬ 
lution of also the differential equations for the structure functions are 

obtained easily [10] 


d 

dL 


< ( v(L)) q >= 


q -1 

£ 


q'- 


{<1 - n)\ 


< D^v q ~ n >. 


(4) 


If the Kramers-Moyal coefficients have the form = d q v q /L (where d q are 
constants) scaling behavior of (1) is guaranteed with £ q = — ]Cn=i (q- n y. dn- For 
a discussion of this finding in terms of multifractality and multiaffinity see [19]. 


4 
























The 1 /L dependence of D ^ indicates that the natural space variable is InL. Thus 
we take in the following as the space variable Aj = ln(Lo/Li). Furthermore we 
normalize the velocity increments to the saturation value of (Too = \J < > for 
length scales larger than the integral length. 

An important simplification is achieved if the fourth Kramers-Moyal coefficient 
is zero, than the infinite sum of the Kramers-Moyal expansion (2) reduces to the 
Fokker-Planck equation of only two terms [18]: 


5a p '»' a > = 




P(v, A). 


(5) 


Now only the drift term and the diffusion term D ^ are determining the 
statistics completely. Form the corresponding Langevin equation, characterizing 
the evolution of a spatially fixed increment, we know that D W describes the 
deterministic evolution whereas D^ describes the noise acting on the cascade. If 
D shows a v-dependency one speaks of multiplicative noise. 

In Fig. 2 we show the evaluated conditional moments M <1> and M <2 ' > as they 
evolve for SX —> 0. For the higher Kramers-Moyal coefficients we found that 
Af( 4 ) < 10~ 2 M^ 2 \ Thus we take the Fokker-Planck equation as the adequate 
description. 




Figure 2 : The conditional moments M^(v,X,5X), (a), and M^ 2 \v, A, SX), (b), 
for A = 1.48 and AA= 0.04 (o), 0.06 and 0.20. The fits to the experimental data 
by polynomial functions (see text) are shown by solid lines. These fits were done 
only with respect to the interval — 2 < v < 2 ). 

To estimate the D ^ and D ^ coefficient properly we have performed the 
following three procedures based on the assumption that we can approximate 
the Kramers-Moyal coefficient by the following polynomials: — 7 (A) and 

£>( 2 ) = a(A) + (5(A) * 1 ^ + /3(A) * vf. At first, we used the same polynomial 
forms for the M-coefficients and estimated the coefficients of the polynomials 
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for S A —> 0. Here problems arise due to the finite value of the Markov-length 
and due to the finite resolution of our detector. Secondly, inserting the polyno¬ 
mials of and D ^ into (4) we can determine directly the coefficients a to <5 
from the measured structure functions (the highest order of the structure func¬ 
tion we used was 6 ). Problems arise due to the noise of the derivative of the 
structure functions, thus only the magnitude of the coefficient could be estimated. 
Thirdly, we have evaluated analogously the corresponding structure functions of 
Peven(v, L) = 1/2 (P(v, L) + P(-v, L)) and P odd (v, L ) = l/2(P(v, L) - P(-v, L)). 
[ 20 ] 

From all these estimations we got a good guess of the values of a(A) to (5(A). 
The best results we obtained for 7 . The worst result was obtained for /?, for 
which we finally used the value n /18 expected from the Kolmogorov picture of 
intermittency [21, 10]. As the final values we got: 

7 (A) = 0.36ea;p 

cx(A) = 0.02e;rp 
(5(A) = 0.013exp(-1.2A) 

0 = 0.017 . ( 6 ) 

Next, let us briefly comment this result. For 7 we find that it has a constant 
value of 1/3 in the inertial range (0 < A < 2). The functional dependence of 7 
is given to extend the A dependence into the dissipation range. The values <5 and 
a are violating scaling behavior. Both, a and <5, are decaying exponentially with 
A. For scales larger than the integral length these two coefficients are large and 
are responsible for the building up of the skewness, because these terms allow the 
change of the sign of the velocity increments within the stochastic process. 

As further test of the validity of our approach we have calculated the evolution 
of the pdfs by a numerical iteration with the Fokkcr-Planck equation. As an 
initial condition the pdf P{vq,Lq) at the integral scale was approximated by an 
empirical function and then inserted into the numerical iteration. We found that 
the evolution of the pdfs depends sensitively on the chosen coefficients. (Finally, 
we changed the numbers in (??) in the range of some percent to obtain a best 
result.) The result of the numerical iteration of the Fokker-Planck equation with 
the above mentioned coefficients are shown in Fig. 1 by solid lines. 

4 Universality 

Let us now discuss how the turbulent cascade will be affected, if the statistics on 
large scales is changed as it may be the case for different driving forces. As an 
extreme case we have chosen a box like form for a large scale pdf (the standard 
deviation and the skewness was adapted to the one of our experimental pdf for Lq 
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(see Fig. 3)). Next, using the above mentioned coefficients ( 6 ), we have iterated 
simultaneously the experimental approximation of P(vo,Lq), see Fig. 1, and the 
box like pdf. As shown in Fig.3 we see that under the iteration the synthetic pdf 
becomes more and more similar to the experimental ones. This convergence has 
been quantified by the y 2 measure (y 2 ex ^(Pi — P 2 ) 2 /Pi )■ We find that \ 2 
decays about exponentially with A. Thus, we see that for a sufficient large cascade 
any large scale pdf will converge against universal small scale pdfs. 



Figure 3: Evolution of a synthetic box-shaped pdf (broken line) under the numeri¬ 
cal integration with the Fokkcr-Planck equation using the coefficients of ( 6 ). From 
top to bottom A = 0, 1.4, 2.2 and 3.0. Presentation as in Fig.l. 


5 Conclusion 

We have presented further evidence that the statistics of fully developed turbulence 
is based on a Markov process in the space variable L. It has been shown how the 
drift coefficient and the diffusion coefficient of a describing Fokker-Planck equation 
can be determined from experimental data. With these coefficients it is possible 
to calculate accurately the pdfs of a cascade. As a further result we have shown 
that situations where the large scale statistics differ, for example due to different 
driving forces, will show the same small scale statistics, as long as the same Fokker- 
Planck equations apply. Therefore, universality of small scale turbulence can best 
be characterized by comparing the drift and diffusion coefficients. 

The next challenging questions are, how do these result depend on the Reynolds 
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number and what happen in anisotropic turbulence. For the ladder point recently 
evidence has been found that for sufficiently small scales again the results of local 
isotropic turbulence hold [22]. We believe that these findings do not only play an 
important role for the complete statistical characterization of turbulence but also 
should have quite practical importance for numerical simulations of flow situations 
with high Reynolds numbers. 
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